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Abstract 


6^ 

jk  recently  developed  solution  for  the  plastic  strain,  e^(x,t) ,  on 
the  crack  line  is  used  in  conjunction  with  a  critical  strain  criterion 
to  construct  curves  for  K^(a)  versus  a,  where  a  is  the  increase  in  crack 

length.  Resistance  curves  have  been  computed  for  various  values  of  the 

Sup  ^ 

critical  plastic  strain.  They  show  a  monotonic  increase  of  K^(a)  with 
increase  in  crack  length,  to  a  constant  steady-state  value. 
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1.  Introduction 


For  small-scale  yielding,  the  condition  for  continued  crack  advance 
may  be  written  as 


K  =  K^a)  , 

where  K  on  the  left-hand  side  is  regarded  as  the  "applied"  K,  and  a  in 


K^Ca)  is  the  increase  in  crack  length.  The  curve  of  ^(a)  versus  a  is 
called  the  resistance  curve  or  R-curve.  Under  plane  stress  conditions, 

K^Ca)  tends  to  be  a  monotonically  increasing  function  of  a,  at  least  for 
many  metals.  Following  the  notation  of  Ref . [l,p.20] ,  we  will  use  K^c  to 
denote  the  value  of  K  for  initiation  of  crack  propagation,  i.e.,  KIC  =  1^(0) . 

In  this  paper  we  consider  a  general  geometry  such  as  an  edge  crack  in 
a  thin  sheet,  and  we  construct  resistance  curves.  These  curves  are  based  on 
a  critical  strain  criterion  for  crack  propagation,  which  stipulates  that 


crack  growth  will  proceed  when  the  plastic  strain  on  the  crack  line  maintains 

Pf 


The 


a  critical  strain  level  €  at  a  distance  xt  ahead  of  the  crack  tip. 

y  f 

critical  strain  criterion  was  proposed  by  McClintock  and  Irwin  [2]  and  it 
has  been  used  frequently  for  the  Mode-Ill  case.  Rice  [3,p.281]  presented  an 
expression  for  the  plastic  shear  strain  on  the  crack  line  for  Mode-Ill  crack 
propagation  in  an  elastic  perfectly-plastic  material.  It  was  also  shown  in 
Ref. [3]  that  the  fracture  criterion  of  critical  plastic  shear  strain  leads  to 
an  integral  equation  for  the  plastic  zone  size  required  for  quasi-static 


crack  extension. 
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The  construction  of  curves  for  K^a)/!^  presented  in  this  paper  for 
plane  stress,  follows  in  general  outline  the  method  of  Ref. [3].  The 
construction  is  based  on  an  expression  for  the  crack-line  strain  which  was 
recently  presented  by  Achenbach  and  Dunayevsky  [4]  and  Achenbach  and  Li  [5] . 
As  shown  in  some  detail  in  Ref. [5],  in  the  plastic  loading  zone  the  stresses 
and  strains  can  be  expanded  in  powers  of  the  distance,  y,  to  the  crack  line. 
Substitution  of  the  expansions  into  the  equilibrium  equations,  the  yield 
condition  and  the  constitutive  equations,  yields  a  system  of  simple  ordinary 
differential  equations  for  the  coefficients  of  the  expansions.  This  system 
is  solvable  if  it  is  assumed  that  the  stress  o  is  uniform  on  the  crack  line 

y 

By  matching  the  relevant  stress  components  and  particle  velocities  to  the 
dominant  terms  of  appropriate  elastic  fields  at  the  elastic-plastic  boundary 

P 

a  complete  solution  was  obtained  for  the  plastic  strain,  e^,  in  the  plane  of 
the  crack.  The  solution  depends  on  position  on  the  crack-line  and  time, 
and  applies  from  the  propagating  crack  tip  up  to  the  moving  elastic-plastic 
boundary . 

It  is  shown  in  this  paper  that  the  critical  strain  criterion  yields 

an  integral  equation  for  x  (a) ,  where  x  (a)  is  the  extent  of  the  plastic 

P  P 

zone  on  the  crack  line  as  a  function  of  the  increase  in  crack  length. 

The  integral  equation  has  been  solved  numerically  for  various  values  of 
Pf 

e  /(£  )r,„»  where  (e  )n_  is  the  crack-line  strain  at  the  elastic-plastic 
y  y  PB  y  PB 

boundary.  A  relation  between  and  x  (a)  subsequently 
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yields  resistance  curves  K^(a). 

The  geometry  is  shown  in  Fig.  1.  The  x^-axis  of  a  stationary  coordinate 
system  is  parallel  to  the  crack  front ,  and  x^  points  in  the  direction  of 
crack  growth.  The  position  of  the  crack  tip  is  defined  by  x^  =  a(t).  A 
moving  coordinate  system,  x,y,z,  is  centered  at  the  crack  tip,  with  its 
axes  parallel  to  the  x^x^  anc*  x3  axes. 
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2 .  Crack-Line  Strain 

For  monotonic  loading,  expressions  for  the  strain  rate  £^(x^,t),  on  the 
crack  line  ahead  of  a  propagating  crack  tip,  are  given  in  Refs.[  4]  and 
[  5  ]  as 


(x  t)  =  -2g.(Q- 

ky  1’  x^-a (t) 


£n 


Xp^ 

x^-a(t) 


B1a(t)+B„x  (t) 

1  2  P 

x1~a(t) 


C..a(t)+C  x  (t) 
_ £_£ _ 

Xp(t) 


[x1~a(t) ] ' 


(2.1) 


where  x^  defines  the  point  of  observation  in  the  stationary  coordinate  system, 
while  a(t)  defines  the  position  of  the  crack  tip  and  x^(t)  is  the  size  of  the  plas¬ 
tic  zone  along  the  crack  line.  Also,  k  is  the  yield  stress  in  shear, 

E  is  Young's  modulus,  and  the  constants  B^,  B^.  and  have  been  derived 
in  Ref . [  4 ]  as 


Bx  =  -|(1+v)[(k  +  5)  +  2(k  +  l)/2]  -  | 

B2  =  |(1+v)[(k  +  5)  +  (k  +  l)/2] 

Cx  =  |(1+v)[-(k  +  5)  +  (k  +  l)/2]  +  | 

C2  =  -|(1+v)[-(k  +  5)  +  2 (k  +  l)/2]  , 
where  v  is  Poisson's  ratio  and  k  is  defined  as  (3-v) / (1+v) . 


(2.2) 

(2.3) 

(2.4) 

(2.5) 


In  the  stationary  coordinate  system,  (2.1)  can  be  integrated  to  yield 


Cy  (Xl’t)  =  ^pb  +  •  (2  -6) 

where  (£y)pg  is  the  elastic  strain  at  the  elastic  -  plastic  boundary: 

(Vpb  *  ■ 


(2.7) 


<*• 
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e^(x^,t)  is  the  plastic  strain 


£y(x 


x>t)  =  /  e  (x1,s)ds 


(2.8) 


The  lower  limit  t  ,  which  is  the  time  that  the  elastic-plastic  boundary 
P 

reaches  position  x^  and  the  plastic  strain  starts  to  accumulate,  follows 
from 


a(t  )  +  x  (t  )  =  x 
P  P  P  1 

By  integration  by  parts  of  the  terms  multiplying  B2  and  C^, 


(2.9) 


e^(x^,t)  can  be  separated  into  two  components: 

t 

e^(x  ,t)  =  eSP[?(x  ,t)]  +  f  epP(x  ,s)ds, 
yx  y  ± 


(2.10) 


where 


and 


|  eJP[C(x1,t)]  =  B2(C-1)  -  \  C2  (^  -  1) 

5(x  ,t)  =  x  (t)/[x  -a(t)]  , 

1  p  1 


^  ^(x^s)  =  a (s)  f[x1-a(s),x  (s)] 


(2.11) 

(2.12) 


(2.13) 


£[Va(s),yS)]  * 


{2£r.C(x1>s)  +  B1~B2C(x1,s)  + 


[((x^S)]3 


[C(x1>S)]2 


(2.14) 


v'  v.v:  .  I'  -S-S  l .  .  .  1 1  y;v;Vy  v'  '  ]  . 
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Here  we  have  used  (2.9)  at  the  lower  limits  of  integration  of  Eq.(2.8). 

jig 

The  lower  limit  t  in  (2.10)  is  defined  as  t  =  t^  for  t^  >  t  ,  where 

tf  is  the  time  that  crack  propagation  starts.  When  t  >  t^ ,  the  lower 

* 

limit  is  defined  as  t  =  t  . 

P 

We  will  now  assume  that  crack  propagation  is  governed  by  the  criterion 


of  a  critical  plastic  strain  at  a  fixed  micro-structural  distance  ahead 
of  the  crack  tip,  i.e.»  at  x^  =  a(t)  +  x^.  The  condition  for  stable 
crack  propagation  then  is 


Pf  sPr 
e  =  e  [x 
y  P 


[x  ( t ) / xf ]  +  J  [a (t)  +  xf,s]  ds. 


(2.15) 


where  is  the  critical  value  of  the  plastic  strain,  and  the  right-hand 

side  follows  from  (2.10).  Equation  (2.15)  is  an  equation  for  x^(t).  For 

the  case  t,.  <  t  ,  the  integral  in  (2.15)  vanishes  when  the  upper  limit  is 
f  P 

taken  as  t  =  tf,  and  the  following  equation  for  x^(t^)  is  obtained 


IT  Xp(tf}  +  [7  C2  '  B2  "  k  Eyf]  XpUf)  ‘  7  Vf  =  ° 

A  convenient  form  of  (2.15)  can  be  obtained  by  replacing  the 


(2.16) 


independent  variable  t  by  crack  length  a.  Thus  we  consider  the  strain  as 
well  as  the  position  of  the  elastic-plastic  boundary  as  functions  of  a. 
After  introducing  the  normalizations 


a  =  a/xf  ,  xp(a)  =  xp[t(a)]/xf  , 


(2.17a,b) 


we  can  write 


v -y  v-*'-  ^iwrwvw  *’.’  *'.’  *‘."ir'  -.’v  t~  ~  •r  ^  -  t-  ^  -  -• 
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e  =  e  [x  (a)]  +  f  f  [a  +  l-n,x  (n)]dn 

y  P  J  *  p 

a 


(2.18) 


where 


c  c 

f[a  +  1  -  r),x  (n)  ]  =  — - —  [2Z.nt,  +  Bn  -  B  ?  +  — - -] 

P  a+l-n  1  2  5 3  C2 


?  =  x  (n)/(a  +  l  -  n) 


(2.19) 


(2.20) 


The  lower  limit  a  in  Eq.(2.18)  is  defined  by 


a  =  0  , 


when 


a  +  1  <  x  (0) 
—  P 


(2.21) 


and 


^  ^  _ 
a  +  x  (a  )  *  a  +  1,  when  x  (0)  <  a  +  1  , 
P  p  — 


(2.22) 


where  x  (0)  =  x  (t-)x-.  Equation  (2.21)  corresponds  to  t  >  t  ,  while 
p  p  t  f  f  p 

(2.22)  corresponds  to  t,  <  t  , 

f  P 

Equation  (2.18)  is  a  Volterra  integral  equation  for  x^(a),  which  can 

be  solved  by  a  step-by-step  procedure.  If  x  (a),  is  known  for 

P 


0  <  a  <  a  <  x  (0)  -  1,  then  for  a  =  a.  +  A  <  x  (0)  -  1,  Eq.(2.18)  yields 

—  —  1  p  _  1  p 

a^+A 

-  J  f[a  +A+l-n>x  (n)]dn  (2.23) 

1  P 


SP ,  Pf 
e  [x  (a.+A)  =  e 

y  p  i  y 


where  A  is  small.  The  integral  in  (2.21)  can  be  approximated  by 
a,+A  a. 


+A+l-n,x  (n)]dn 
p 


.  r  ^ 


+A+l-n.x  (n)]dn  +  f(a  +A+l-n,x  (a)]A  (2.24) 

P  1  P 


o  r  o 

The  integral  on  the  right-hand-side  of  (2.24)  is  known.  Substitution  of 

(2.24)  in  (2.23)  yields  a  cubic  equation  for  x  (a.+A),  which  can  be  solved. 

P  1 

This  value  for  x^(a^+A)  is  used  as  the  starting  point  for  an  iteration  procedure 

whereby  subsequent  values  of  x  (a.+A)  are  substituted  in  the  integral  in 

P  1 


Eq.(2.21)  to  yield  improved  values  of  x  (a.+A),  until  a  desired  accuracy 

P  1 

has  been  achieved.  For  the  first  step  in  this  procedure  we  use 


x  (A)  s  x  (0)  +  x’(0)A,  where  x  (0)  follows  from  (2.16),  and 
P  P  P  P 

x^(0)  =  dXp(a)/da  at  a  =  0,  can  be  obtained  from  (2.18). 

For  the  case  when  x^(0)  -  1  <_  a,  defined  by  (2.22),  the  lower  limit  of 
the  integral  is  a  function  of  the  upper  limit,  a,  and  of  the  unknown  function 

x^(a).  Again,  if  x^ (a)  is  known  for  a  a^,  then  for  a  =  a^  +  A ,  Eq .  (2.22) 

*  -  *  —  _  *  *  _ 
yields  a  +x(a)=a1+A+l,  where  x  (a  )  is  known  since  a  <  a.  Hence 

pi  p 

*  *  — 

we  can  solve  a  =  a  (a^,A).  Using  an  analogous  method  as  before,  Eq.(2.18) 

can  subsequently  be  solved  for  x  (a  +A) . 

P 

In  the  actual  computations  the  following  values  of  the  relevant  parameters 
have  been  considered: 


ey  E  eyf/(ey}PB  =  2’  6,and  10  (2.25) 

The  Poisson's  ratio  was  taken  as  v  =  0.3.  For  a  =  0,  the  value  of  x  (0) 

P 

follows  directly  by  computation  of  the  real  root  of  Eq.(2.16),  since 
Xp(0)  =  Xp(t^)/X£.  The  numerical  results  show  a  subsequent  increase  of 
Xp(a)  with  a,  where  x^(a)  approaches  an  asymptotic  value  for  large  a. 

For  quasi-static  steady-state  crack  extension  the  following  relation 
was  derived  by  Achenbach  and  Li  [5]: 

£y (x)  {[*n(^-)]"  -  B^n(J-)  -i  [(^-)3  -  l]}  (2.26) 

P  P  P 

The  critical  strain  criterion  then  yields  the  relation 

f  ~  [£n(x  )]2  -  B  tax  +  [0/x  )3  -1]  =  0 

K  y  p  1  p  3  p 


(2.27) 


Equation  (2.27)  can  be  solved  for  to  yield  a  result  which  is  independent 

of  crack-tip  speed  and  loading  state.  This  result  is  just  the  asymptotic 

value  of  x  (a)  as  a  increases. 

P 

A  convenient  form  of  the  resistance  curve  shows  the  ratio  K^/K^  versus 
a  dimensionless  crack  length.  Here  is  the  present  stress  intensity 
factor,  and  K^c  is  the  value  of  the  stress  intensity  factor  which  is 
required  to  satisfy  the  fracture  criterion  for  a  stationary  crack.  In 
Ref. [5,  Eq.57]  an  expression  was  presented  which  relates  x^(t)  to  K^/k. 

In  the  present  notation  this  expression  states 


xp(a) 


2/2  1  1  ,  ,  , 2 

[Ki(a)/kl 


It  follows  that 


(2.28) 


K1(i)/KIC  =  KR(a)/KIC  =  [x  (5)/x  (O)]1*  (2.29) 

The  quantity  x^  which  enters  in  a,  can  be  eliminated  by  equating  the  result 

for  x  (0)  obtained  from  (2.16)  to  x  (0)  as  obtained  from  (2.28).  By  using 
P  P 

the  resulting  expression  for  x^  in  the  definition  of  a,  as  given  by  (2.17a), 
we  f ind 


a 


9tt 


a  x  (0) 

_ E _ 


2/2  [KIC/k]2 


(2.30) 


For  the  three  values  of  given  by  Eq.(2.25),  curves  for  K^(a 

versus  a  are  shown  in  Fig.  2,  It  is  noted  that  the  resistance  curves  show 
a  monotonic  increase  with  increase  in  crack  length,  to  a  stable  phase  of 
crack  propagation  where  K^(a)/Kp^  assumes  the  steady-state  value. 
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